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Abstract. - Fluids adsorbed in 3D wedges are shown to exhibit two types of continuous 
interfacial unbinding corresponding to critical and tricritical filling respectively. Analytic solu- 
tion of an effective interfacial model based on the transfer-matrix formalism allows us to obtain 
the asymptotic probability distribution functions for the interfacial height when criticality and 
tricriticality are approached. Generalised random walk arguments show that, for systems with 
short-ranged forces, the critical singularities at these transitions are related to 2D complete 
and critical wetting with random bond disorder respectively. 



Recent studies of filling transitions for fluids in 3D wedges [1,2] have revealed the much 
stronger influence of interfacial fluctuations compared with wetting at flat and rough sub- 
strates [3-7]. Encouragingly effective Hamiltonian predictions for the critical exponents at 
continuous (critical) wedge filling with short-ranged forces have been confirmed in large scale 
Ising model simulation studies [8] . Similar experimental verification of the predicted geometry- 
dominated adsorption isotherms at complete wedge filling [9] raise hopes that the filling tran- 
sition itself and related fluctuation effects will be observable in the laboratory. Here we further 
develop the fluctuation theory of 3D filling and show that there is a rather deep and previously 
unrecognized connection with the theory of wetting in 2D systems with random bond (RB) 
disorder. Our findings are based on analytical solution of an effective model of 3D wedge 
filling and also generalised random- walk arguments [10]. First we show that there are actually 
two types of continuous filling behaviour corresponding to critical and tricritical transitions 
respectively with the latter having stronger fluctuation effects. The phase diagram for these 
transitions together with the classification of fluctuation regimes and the allowed values of 
critical exponents resemble very closely those predicted for 2D wetting. More precisely 3D 
critical filling is related to 2D complete wetting whilst 3D tricritical filling is related to 2D crit- 
ical wetting. Remarkably the particular value of the 3D wedge wandering exponent for pure 
systems (thermal disorder) implies that criticality at tricritical and critical filling is related to 
predictions for 2D critical and complete wetting with RB disorder. 

Consider the interface between a bulk vapour at temperature T and saturation pressure 
with a 3D wedge characterised by a tilt angle a (see Fig. ^1 . Macroscopic arguments dictate 
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that the wedge is partially filled by liquid if the contact angle 6 > a and completely filled 
if 9 < a [11] and is fully supported by both interfacial Hamiltonian [12] and exact Ising 
studies [13]. The filling transition refers to the change from microscopic to macroscopic liquid 
adsorption as T — ► Tf, at which 0(Tf) — a, and may be first-order or continuous (critical 
filling). Both of these transitions can be viewed as the unbinding of the liquid- vapour interface 
from the wedge bottom. Characteristic length scales are the mean interfacial height above 
the wedge bottom lyy, the roughness £j_ and the longitudinal correlation length £ y , measuring 
fluctuations along the wedge (see Fig.^). The relevant scaling fields at critical filling are 8 — a 
and the bulk ordering field (partial pressure) h. In our discussion of filling we shall work 
exclusively at bulk coexistence (h = 0) since it is here that the connection with RB wetting 
emerges. However calculations away from coexistence are not in any way problematic. At 
coexistence we define critical exponents by lw ~ {6 — a)^ l3w and £ y ~ (9 — ct)~ Vy . The 
roughness can be related to £ y by the scaling relationship: 

£t ~ 4 w (i) 

where Qw is the wedge wandering exponent. 

The liquid-vapour interface across the wedge is aproximately flat and soft-mode fluctua- 
tions arise from local translations in the height of the filled region along the wedge [1] . The 
pseudo-one-dimensional nature of these means that Cw is greater than the wandering expo- 
nent defined for the asymptotically fiat free interface £3 in a 3D system. For systems with 
sufficiently short-ranged forces critical filling is fluctuation-dominated (i.e. £j_ ~ lw) an d 
dimensional reduction arguments lead to the identification [2] : 

1 + (2 - (3 

with Qi the 2D free interface wandering exponent. Thus for pure (thermal) systems for which 
Cd = (3 — d)/2 for d < 3, the 3D wedge wandering exponent C,w = 1/3. In the fluctuation- 
dominated regime the value of Cw determines the other exponents, in particular 

(3) 



2C 



w 



with Pw — v ±- — Cw v y Thus for thermal forces v y = 3/4, a value verified in Ising simulation 
studies [8]. 

There is, however, another example of continuous filling with even larger fluctuations 
characterized by different scaling fields and critical exponents. This corresponds to tricritical 
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filling and there are two ways the transition may be induced. The first mechanism occurs for 
wedges made from walls which themselves exhibit weakly first-order wetting transitions [1]. 
The filling transition in these systems will be first-order if a < a* and critical for a > 
a*, where a* is the tricritical wedge angle. A second mechanism exists, however, which is 
more practicable in simulation studies. Imagine a wedge made from a homogeneous chemical 
material which exhibits critical filling. Now micropattern a stripe along the wedge bottom (see 
Fig-D), so as to weaken the local wall- fluid intermolecular potential and therefore strengthen 
the interfacial binding potential, since locally liquid adsorption will not be favoured and 
consequently the vapour-liquid interface is more likely to be localized near the wedge bottom. 
This situation can be easily engineered in Ising model studies by modifying the strength of 
the surface field near the wedge bottom and indeed has been done on planar substrates in 
the laboratory. With this modification it may be possible to bind the interface to the wedge 
bottom even at the filling boundary 8 = a and at bulk coexistence. A continuous tricritical 
transition may then be induced as the strength of the modified wedge potential approaches a 
tricritical value. 

For thermal systems, both critical and tricritical transitions can be modelled by the wedge 
Hamiltonian [1] 

H w [l } = Jdyf^^y + VwQo)} (4) 

where lo(y) is the local height of the interface at position y along the wedge bottom and 
£ is the liquid-vapor surface tension. This expression arises from the identification of the 
interfacial breather modes which translate the interface up and down the sides of the confining 
geometry, as the only relevant fluctuations for the critical filling transition, and corresponds 
to the excess free-energy contribution of a breather-mode configuration with respect to the 
planar case obtained from the usual capillary wave model [1] . Note that the effective bending 
term resisting fluctuations along the wedge is proportional to the local interfacial height since 
the interfacial cross-section has a length 2lo/a. At bulk coexistence, the effective binding 
potential Vw{Iq), up to unimportant additive constants, is given by [1]: 

VwQo) « ^° 2 - a2 ^ + AV w (l Q ) (5) 
a 

where AVw(lo) has a hard wall repulsion for Iq < and a long-ranged tail which decays 
when Iq — > oo. For short-ranged forces this can be modelled as a contact-like potential with 
strength u (i.e. AVw ~ —uS(lo), where 5(x) is Dirac's delta function). In general, there will be 
a tricritical value u c such that for u < u c the interface unbinds from the wedge bottom when 
8 — > a whilst for u > u c , the interface remains bound to the wall in the same limit. A section 
of the phase diagram (at h = 0) for this is presented in Fig. Ota) and shows two continuous 
filling transitions. Critical filling corresponds to 8 — > a for u < u c (route (iii)) and has the 
critical exponents described above. Tricritical filling corresponds to any thermodynamic path 
for which u — > u c and — > a (routes (i) and (ii), for example). Along the path 8 = a the 
relevant length scales diverge as: 

l w ~ (u - u c )-^' , £y ~ (u - u c )~ v 'y (6) 

and £x ~ t,y w ■> with the tricritical wandering exponent. These expressions define new 
critical exponents which are distinct from those at critical filling. Again for short-ranged forces 
we anticipate that the transition is fluctuation-dominated with lw ~ and (3^ = Cw v y- 
More generally, in the vicinity of the tricritical point (and at h = 0) we anticipate scaling 
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Fig. 2 - Phase diagrams for (a) filling and (b) wetting transitions in terms of the ordering fields 9 — a 
and h and the contact-like potential strengths u and w. The thick and dashed lines in both diagrams 
correspond to continuous and first-order boundaries between bound and unbound interfacial states, 
respectively. The arrows show representative paths along which continuous unbinding occur. The 
filled circles represent the tricritical filling and critical wetting points, respectively. 

e.g. £, y <~ \u — u c \~ v y A [(9 — a)\u — u c |~ A *] with the gap exponent A*. Thus along route (ii) 
S v ~(0-a)-W. 

To continue we evaluate these exponents for thermal systems by explicit transfer-matrix 
analysis before recasting them more generally in terms of the wandering exponent (w In the 
continuum limit the partition function is defined as a path integral [14]: 



z[h,i a ,Y] = J vi cM-n w [io]) 



(7) 



where Y is the wedge length, l a = ?o(0) and If, — lo(Y) are the endpoint heights and we have 
set ksT = 1 for convenience. Due to the presence of a position-dependent stiffness some 
care must be taken with the definition of the path integral. This turns out to be of crucial 
importance for the the evaluation of the exponents at tricritical (but not critical) filling. This 
problem was already pointed out in Ref. [15] and is related to the well-known ordering problem 
in the quantization of classical Hamiltonians with position-dependent masses. Similar issues 
also arise in solid state physics [16]. Borrowing from the methods used to overcome these 
difficulties we use the following definition 



N 

Z[l b ,l a ,Y}= lim / dl 1 ...dl N - 1 T[K(l j ,l j _ 1 ,Y/N) 



where Iq = l a and Zjv = h, and K(l, I', y) is defined as: 



K(l,l',y) 



exp 



[l-l'f -yV w (l) 



airy \ ay 
In the continuum limit the partition function becomes 

Z{l b ,l a ,Y)=Y J ^n{lbWn{Qe- E ^ 



where the complete orthonormal set of functions satisfies 

3a 
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\ip = Eip 



(8) 
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In the thermodynamic limit Y — > oo we obtain the probability distribution function (PDF) 
for the midpoint interfacial height PwQ-o) — iV'o^o)! 2 an d the longitudinal correlation length 
£y = l/(Ei — Eq). We have obtained the analytical solution to the transfer-matrix operator 
for short-ranged forces determining the crossover from tricritical to critical filling. We present 
without proof some of our findings. A more complete and detailed description of our results 
will be published elsewhere. 

Along route (i) we find that there is only one bound solution to Eq. for u > u c with 
Eq cx (u — u c ) 3 and associated PDF 

P w (l Q )Kl Ai 2 [c(u-u c )l } (12) 

with c an unimportant constant while Ai(x) is the Airy function. Thus Vw ~ £_L oc (u — u c ) _1 
and £ y cx (it — u c )~ 3 identifying (3^ = 1, v* = 3 and confirming that the tricritical wandering 
exponent coincides with the critical wandering exponent = Qw = 1/3. A numerical plot 
of the PDF is shown in Fig. [31 On the other hand, the scaling of the PDF along the route (ii) 
is given by 



, '-'< h 2/;; 

Pw{h) oc Z exp — - T 

4e 4e 



where £ e = 5T 1/2 [(6</a) 2 - l]" 1 / 4 , s = e 2 /4 - 1/2 with e = -Y,E Q/a « 1.086 and H s (z) 
is the Hermite function [17]. The value of e is obtained by imposing appropriate boundary 
conditions at 1$ = 0. Thus along this route lw ~ £j_ oc (6 — a) -1 / 4 similar to critical filling. 
From analysis of the spectrum it is also possible to show that £ y cx (6 — a)~ 3 / 4 . As anticipated, 
in the vicinity of the tricritical point the divergent length scales show scaling with tricritical 
gap exponent A* = 4. Whilst the exponents for critical filling are already known the exact 
scaling form for the PDF has not been given previously. For thermodynamic paths (iii) far 
from the tricritical point we have found that the scaling of the PDF is of the form shown in 
Eq. Q13[) but with e w 1.639. We remark that the PDFs at critical and tricritical filling have 
distinct short-distance expansions when Iq/Iw ^ and our results (cubic and linear powers 
respectively) are consistent with exact thermodynamic requirements [2]. 

We can now place the results for tricritical filling in a more general context. For short- 
ranged forces the tricritical transition belongs to a general class of strong-fluctuation regime 
interfacial unbinding since the transition occurs at a finite value of the binding potential u. 
The critical singularities at such transitions can be very elegantly modelled using random- 
walk methods which previously have been succesfully used to understand 2D critical wetting 
(at flat walls) [10]. By modelling the interface as a sequence of bound and unbound regions 
(the so-called bead-necklace picture) , the critical exponent for the correlation length along the 
interface at the strong-fluctuation transition can be related to the appropriate wandering expo- 
nent. Details will be published elsewhere, so we just mention our main results. Following this 
argument for the tricritical filling transition and making allowance for the position-dependent 
stiffness, we find 

v* = 1 , A* = 2 ~ 2Cw (14) 

from which all other critical exponents follow. For thermal systems, Qw = 1/3, implying 
v* — 3 and A* = 4, in agreement with our explicit calculations. Eqs. (|12|) . (I13fl and l|14|) 
are the main results of our paper and together with Eq. @ completely determine the critical 
singularities at fluctuation-dominated 3D filling occuring at bulk coexistence. 

A remarkable connection with the theory of 2D complete and critical wetting is now 
apparent. These transitions correspond to the continuous unbinding of an interface from a 
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Fig. 3 - Scaled interfacial height PDF Pw(c(u — u c )Iq) as a function of the scaled wedge midpoint 
interfacial height c(u — u c )lo for 9 = a. Inset: Plot of the scaled PDF along routes (ii) and (iii) in 
Fig-IHa), i-e. for e w 1.086 (continuous line) and e « 1.639 (dot-dashed line), respectively. 



planar wall and can be modelled by the interfacial Hamiltonian 



nil] = / dx 



(15) 



where l(x) is the interfacial height at a position x along the wall, £ is the 2D stiffness and 
W(l) is the binding potential. In general W(l) — hi + AW (I) where h is proportional to the 
bulk ordering field (partial pressure) and for short-ranged forces AW(l) can be modelled as 
a contact potential of strength w (i.e. AW (I) = —wd(l)). Disorder arising from bulk random 
bonds can also be allowed for by including a stochastic term in W(l). The phase diagram 
is shown in Fig. H£b), and shows two continuous transitions referred to as complete wetting 
(path (hi)) and critical wetting (e.g. path (i) and (ii)), at which the mean interfacial height 
l n , roughness £j_ and transverse correlation length £|| all diverge. For short-ranged forces 
the transitions are fluctuation-dominated and i w ~ £j_ ~ £jj 2 . More specifically, for complete 

wetting we write £|| ~ h~ v W whilst for critical wetting £|| ~ (w — w c ) _l/ ii with associated gap 
exponent A off coexistence. The values of these exponents can be expressed explicitly in terms 
of the wandering exponent (2 as [18] 

There is therefore a qualitative and quantitative connection between critical singularities at 
3D filling and 2D wetting. The phase diagrams are equivalent (see Fig.^J with the field 8 — a 
playing the role of the ordering field h at wetting. Further, writing the 2D wetting exponents 
in terms of (2, e.g. = ^"(i^), etc., we have from Eqs. © and (|14f> the dimensional 
reduction relations for the critical exponents 

v y = v^{2C, w ) , v* y = is\\(2Cw) , A* = A(2( w ) (17) 

Most remarkably for thermal forces the numerical value 2£jy = 2/3 means that 3D filling is 
related to 2D wetting with RB disorder, since C2 = 2/3 for this case [19]. 
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A full account of our study, including the analysis of the presence of long-ranged forces, 
will be presented elsewhere. 

In summary we have identified a second example of continuous filling transition corre- 
sponding to tricritical behaviour. For systems with short-ranged forces and thermal disorder 
only we have exactly found the critical singularities and associated probability distribution 
function for the interfacial height. A random- walk analysis reveals a remarkable connection 
between the critical exponents for thermal 3D filling and random-bond 2D wetting for systems 
with short-ranged forces. These predictions may certainly be tested in Ising model simulation 
studies and would be a stringent test of the theory of 3D wedge filling. Experimental studies 
of tricritical (and critical) filling similar to those already performed for complete filling would 
be very welcome. 

* * * 

J.M.R.-E. acknowledges financial support from the European Commission under Contract 
MEIF-CT-2003-501042. 
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